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Coherent control of collective spontaneous emission in an extended atomic ensemble resonantly 
interacting with single-photon wave packets is analyzed. A scheme for coherent manipulation of 
collective atomic states is developed such that superradiant states of the atomic system can be 
converted into subradiant ones and vice versa. Possible applications of such a scheme for optical 
quantum state storage and single-photon wave packet shaping are discussed. It is shown that also 
in the absence of inhomogeneous broadening of the resonant line, single-photon wave packets with 
arbitrary pulse shape may be recorded as a subradiant state and reconstructed even although the 
duration of the wave packets is larger than the superradiant life-time. Specifically the applicability 
for storing time-bin qubits, which are used in quantum cryptography is analyzed. 
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I. INTRODUCTION 



Motivated by the development of modern quantum in- 
formation science there is significant interest in inves- 
tigating the interaction of nonclassical states of light 
with atomic ensembles. One active research area is the 
development of optical quantum memories. Such de- 
vices, which can store and reconstruct quantum states of 
light, form a basic ingredient for optical quantum com- 
puters, dealing either with discrete quantum variables 
(qubits) represented by single-photon two-mode states 
of the electromagnetic field or with continuous quantum 
variables represented bv the quadrature amplitudes of the 
field (see reviews P, Q respectively) . There are experi- 
mental demonstrations based on electromagnetically in- 
duced transparency Q for the discrete variable approach 
and off-resonant interaction of light with spin polarized 
atomic ensembles ^ for the continuous variable case. 
The photon echo based quantum memory proposal has 
also attracted considerable attention 0, 0, la 13 • An- 
other cause of the interest is the possibility of atomic 
ensemble state manipulation using atom-photon corre- 
lation (see review The creation of a robust en- 
tanglement of atomic ensembles using Raman scatter- 
ing may be especially important for applications involv- 
ing long-distance quantum communication 'Tl'|. Finally, 
much effort has been directed toward implementation of 
controllable sources of nonclassical states of light such as 
photon-number or Fock states, which represent an essen- 
tial resource for the practical implementation of many 
ideas from quantum information. These sources can also 
be constructed using quantum storage techniques in op- 
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tically dense atomic media |l2j . 

The motivation of all these investigations is that pho- 
tons can interact much more strongly with ensembles 
containing a large number of atoms than with individ- 
ual atoms. In addition there can be a further collective 
enhancement taking place because the absorption, emis- 
sion or coherent scattering processes employed in sev- 
eral of these implementation schemes involve collective 
atomic modes. On the other hand, whenever the num- 
ber of atoms is large, the contribution from the collec- 
tive spontaneous emission jlj, (l^ to the evolution of 
an atomic-field system can also be considerable and may 
need to be taken into account. Therefore the investiga- 
tion of this phenomenon in various situations when light 
interacts with atomic ensembles is of interest. 

In this paper we investigate the possibilities of coher- 
ent control of collective spontaneous emission in an ex- 
tended atomic ensemble interacting with a single-photon 
wave packet and vacuum modes of the electromagnetic 
field. Spontaneous emission is normally considered as a 
noise, which should be inhibited. Various schemes for in- 
hibiting or reducing spontaneous emission noise for the 
single atom case have been proposed and investigated. 
For example in ,16] excitation by 27r-pulses on an auxil- 
iary transition is used to change the phase of the excited 
state wave function and in this way inhibit the sponta- 
neous decay. In an atomic ensemble confined to a volume 
with dimensions small compared to the wavelength of the 
emitted radiation the problem is reduced to the creation 
of antisymmetric subradiant states 0, 0, 0|, which 
are also interesting as they may be a basis for forming a 
decoherence free subspace (19| for the excited ensemble. 
In extended atomic systems, which is a much more fre- 
quent experimental situation, the collective spontaneous 
emission is characterized by a sharp directedness in space 
and occurs only in a small fraction of all radiation modes. 
Therefore it is possible to implement coherent control of 
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collective atomic states using coherent excitation through 
other non-collective modes. Such a manipulation of col- 
lective states may be used for quantum memories in the 
regime of optical subradiance inhibiting the normal spon- 
taneous decay, an idea which was briefly outlined in ,2'ol| . 

The paper is organized as follows. In Sec. II, we intro- 
duce a physical model of interaction between an extended 
many atom system and a single-photon wave packet and 
present an analytical description of the superradiant for- 
ward scattering. In Sec. Ill, we present a scheme for co- 
herent manipulation of collective atomic states that can 
transform superradiant states into subradiant ones and 
vice versa. In Sec. IV, we discuss possible applications of 
such a scheme for the creation of optical quantum mem- 
ory devices and controlled sources of non-classical light 
states. In Sec. V, we propose some specific experimen- 
tal implementations that may be used to observe these 
effects and some possible methods for subradiant state 
preparation. 



II. PHYSICAL MODEL 

Consider a system of ^ 1 identical two-level atoms, 
with positions (j = 1, . . . , N) and resonance frequency 
ujQ, interacting among themselves and with the external 
world only through the electromagnetic field. Let us de- 
note the ground and excited states of jth atom by \0j) 
and The Hamiltonian of the system, in the interac- 
tion picture and rotating-wave approximation, reads 

^ = E ^5k,s^^ak,. e''"'-^ e'(-o-)t + h.c. (1) 



Here 



h \2enV 



1/2 



(d • ek,s) 



(2) 



is the atom-field coupling constant, bj = |Oj)(lj| is the 
atomic transition operator, ak,s is the photon annihila- 
tion operator in the radiation field mode with the fre- 
quency Lo — kc and polarization unit vector ek,s (s — 
1,2), V is the quantization volume of the radiation field 
(we take V much larger than the volume of the atomic 
system), d is the dipole moment of the atomic transition. 
For the sake of simplicity we assume that the vectors £k,s 
and d are real. 

To describe the interaction of the field with the atoms 
we use the approach |2i,.22j. First, it is convenient to 
assume that 

(i) the atomic system has a shape of a parallelepiped 
with the dimensions L^, (a — x, y, z) and the atoms are 
placed in a regular cubic lattice, so that N — N^NyN^ 
and d — LaN~^ is the interatomic distance, 

(ii) the dimensions are much larger, but the inter- 
atomic distance d is much smaller, than the wavelength 
A = 27rcW(^^ of the atomic transition. 



(iii) the center of the system lies in the origin of the ref- 
erence frame. Then we can define the following collective 
atomic operators: 



N 



(3) 



where — 2TmaL^^, Ua — 0, 1, . . . ,Na — 1. Hamilto- 
nian ^ can now be expressed in terms of the collective 
operators obtaining 



q,k,s 



*(q-k) e'('^«~'^)* + H.: 



(4) 



where 0(x) — J^j exp(ix • r^) is the diffraction func- 
tion. The presence of this function in the Hamiltonian 
underlines the fact that each atomic mode q is coupled 
only to modes k lying in a diffraction angle around q (21| . 

We are interested in the interaction of the atomic sys- 
tem with a single-photon wave packet. Therefore we have 
the following general form of the state of the system at 
the initial time t = 



|V'(0))=5]/k,.(0)|0)|lk,.)+5]cq(0)|le 
k,s q 

with normalization condition 



(5) 



(6) 



where |0) — |0i, O2, . . . , Oat) is the ground state of the 
atomic system, |vac) is the vacuum state of the radiation 



field. Ilk . 



Ill J vac) and llr 



7V-i/2i?t|o). 



Now, we assume that 

(iv) a single-photon wave packet propagates in the z- 
direction, the excitation volume may be approximated 
by a cylinder with the cross section S and the length L^, 
and the wave front of the packet is planar inside the ex- 
citation volume, 

(v) the Fresnel number of the excitation volume F = 
S{LzX)~^ « 1, so that the collective interaction of the 
whole ensemble of atoms with the quantum field takes 
place only for longitudinal modes q, 

(vi) the duration of the single-photon wave packet as 
well as that of collective spontaneous emission are much 
greater than cooperative time Tc which in turn is 
much greater than the correlation time of the vacuum 
reservoir te = LzC"^ (Born-Markov condition). Here 



1/2 



(7) 



where Ti is the inverse of the Einstein A coefficient. With 
these assumptions the interaction of the atomic system 
with the electromagnetic field may be described in a one- 
mode approximation with respect to the atomic system 
and a one-dimensional approximation with respect to the 
field. Substituting Eqs. Q and jSJ in the Schrodinger 
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equation and omitting the index q for the single atomic 
mode we obtain 



dt 

dcjt) 
dt 



iV,^(q-k)c(t)e-^("''-")*, (8) 



= -iVN 



(9) 



Now it is convenient to define a single-photon dimcn- 
sionless photon density at the origin of the reference 
frame (r = 0). For the incoming wave packet we have 



1/2 



(10) 



and for the emitted radiation we have the analogous 
equation with Fin{t) and /k,s(0) replaced by F{t) and 
/k.s(0; respectively. Then the solution of Eqs. ((Sj) and 
may be written as (see Appendix 

c{t) = c(0)e-*/2-« 

/ dTFi„(t-T)e--/2-«, (11) 

/tttt 

(12) 



Fit) = F,^it) + J^c{t), 
V '^R 

where tr is the superradiant life-time: 



1 

TR 



iV/i 



1 



1 



47ren ihc^ 



(13) 



and /i = 3A^(87rS')~^ is a geometrical factor, which is 
approximately equal to the ratio of the diffraction solid 
angle of collective emission to Att. If we consider the 
case when c(0) = and substitute Eq. ((TT)l into ((T^ . 
we obtain a solution for superradiant resonant forward 
scattering of photons, which is well known in the the- 
ory of propagation of coherent pulses through a resonant 
medium and especially in the coherent-path model of nu- 
clear resonant scattering of 7 quanta |24| : 



1 

F{t)=F4t) / drFi„(t 

TR Jo 



-t/2tr 



(14) 



Under the approximations mentioned earlier, Eq. 1)14(1 de- 
scribes the simplest regime of superradiant forward scat- 
tering, when the incoming field is scattered only once. 
Considering resonant propagation, e.g. of small area 
pulses as in '25^ , this corresponds to samples with aL not 
far above 1. It may be noted that the collective superra- 
diant atomic state |lq) is known in the theory of nuclear 
resonant scattering as a nuclear exciton (see, for exam- 
ple, H^li^l and references therein). The principal point 
here is the indistinguishability of atoms in a macroscopic 
sample with respect to photon absorption and emission 



into the longitudinal modes of the electromagnetic field, 
which leads to enhancement of atomic-field interaction 
and directedness of collective spontaneous emission j28l |. 
Since the photon propagation time te is much less than 
the cooperative time Tc — y/TRjE, which in turn is much 
less than the duration of incoming and scattered photons, 
propagation effects may be neglected. In fact, from the 
view point of propagating fields, the sample looks like a 
scattering center, which is characterized by an enhanced 
cross-section of the resonant transition with respect to 
longitudinal modes. But the interaction of the atomic 
system with transverse modes of electromagnetic field is 
not collective and these modes do not see an enhanced 
cross section. This is of key importance for the prepara- 
tion of subradiant states in an extended atomic system, 
which will be considered in the next section. 



III. SUPERRADIANT AND SUBRADIANT 
STATES OF AN EXTENDED ATOMIC 
ENSEMBLE 

In his basic paper Dicke considered the two regimes 
of collective spontaneous emission of photons: superra- 
diance and subradiance, which result from the construc- 
tive and destructive interference of atomic states, respec- 
tively. In the first case, a system of inverted atoms under- 
goes the spontaneous transition to the ground state for a 
time inversely proportional to the number of atoms, while 
in the second case the rate of collective spontaneous emis- 
sion, on the contrary, decreases compared to the rate of 
spontaneous emission of single atoms. In the ideal case 
of a localized system confined to a volume with dimen- 
sions small compared to the wavelength of the emitted 
radiation, the rate of collective spontaneous emission of 
photons is equal to zero if the atoms are in an antisym- 
metric collective state. If an atomic system is extended, 
the rate of collective spontaneous emission can be sup- 
pressed only for a few collective modes. In particular, 
one-mode subradiance can be observed in samples hav- 
ing a shape of a cylinder with proportions defined by 
the Fresnel number F 1. It is well known that such 
a one-mode model of collective spontaneous emission is 
equivalent to the Dicke model (2^ , except that the rate of 
spontaneous emission decreases by the geometrical factor 
fi. If we redefine atomic states multiplying them by 
phase a factor exp(zq • Vj) and denote a collective one- 
mode atomic state |nq) corresponding to the excitation 
of n atoms as In), then 



-^i?t|0)==-^^|0,...,l„...,0), (15) 

V ^ i 

:(i?t)2|0) 



|2) = 



^2N{N~l) 



N{N -1) 



-^|0,...,1.,...,1„...,0),(16) 
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The rate of photon emission in this case is {N — 2)fj,T^^, 
i.e., in fact essentially the same as the superradiant spon- 
taneous decay from the state |1) at large N. To convert 
the state in Eq. 1)21(1 into a subradiant state, it is nec- 
essary to irradiate half of the atoms, whose phase was 
changed after absorption of the first photon, and half 
atoms, whose phase was not changed, by additional 2tt- 
pulses. As a result, the atomic system will be divided 
into four parts, which we denote A, B, C, and D, and 
the state of the system will take the form 



FIG. 1: (Color online) Scheme of the working atomic levels, 
which are used for preparing subradiant collective states. 



etc., where R — '^jLi^j- The rate of collective spon- 
taneous emission of photons fiTi~^\{n — l\R\n)\'^ is in- 
creased during excitation of the medium (for n < 
The spontaneous emission of a photon from the state |1) 
occurs with the rate N^T^^ (see Eq. ()13(l 'l. from the state 
|2) with the rate 2{N - l)Airf \ etc. 

We assume now that atoms have the additional level 
\2j) and the transition frequency Wi2 differs from luqi 
(Fig.nj. Then, irradiation by a short coherent 2tt pulse at 
the frequency lui2 can almost instantly change the phase 
of the states thereby changing the phase of the col- 
lective state |1) to the opposite one, while the phase of 
the state |2) remains unchanged. 

Let us divide a sample into two parts, A and i?, each 
containing atoms. It is convenient to represent super- 
radiant states H15() and H16|l in the form 



\U+b) 
|2a+b> 



V2 



(|1a,0b) + |0A,li3)), 



(17) 



1 



2x/lV^ 
+V2N\1a,Ib 



[VN^i\2A,0B) + \0A,2B)) 



(18) 



Upon irradiation of atoms in, for example, part B, by a 
coherent 27r-pulse the atomic system passes from states 
(|17|l and (|18|l to the antisymmetric single-mode subradi- 
ant states 



Il4- 



1 



A-B 



V2 



\2a- 



(|U,0s)-|0a,1s)), 
1 



(19) 



B) = 



2y/N^ 

-V2N\1aAb 



[VN^{\2a,0b) + |0a,2b)) 



(20) 



The rate of photon emission from state H19|) is identically 
zero, while that from state is 2{N - l^-fiT^^, i.e., 
goes to zero for large N. 

Now assume that the atomic system in state H19II ab- 
sorbs a photon and passes to the state 



I'^a-b+c-d) = —^H^a-b^Oc+d) 



\0a+b,2c-d))- 

(22) 



The rate of collective spontaneous emission from this 
state is equal to 4(A^ - 2)" V^T^ 

Therefore, by irradiating the system of atoms by 27r 
pulses through auxiliary transition after each absorption 
of a photon, we can transfer it into the next higher ex- 
cited subradiant state, where the spontaneous transition 
from this state to the previous subradiant state will be 
forbidden. When transfering the atomic system into sub- 
radiant states repeatedly the action of the 2t: pulses can 
be described using Hadamard matrixes 



1 1 

1 -1 



H. 



2k 



(23) 



where k = 2,4,8, etc. The scheme for the preparation 
of different subradiant states can be seen directly from 
the matrix where each column correspond to one specific 
spatial part and, with the exception of the first row, each 
row shows the phase of all the spatial parts for a spe- 
cific subradiant state. Transformation from one row to 



another is accomplished by phase shift \lj 



II,) for 



half of the atoms. The size 2k is equal to the number 
of the spatial parts and the total number of orthogonal 
subradiant states which can be prepare in such a way is 
2A:- 1. 



IV. QUANTUM MEMORY APPLICATIONS 
AND SINGLE-PHOTON PULSE SHAPING 

Now, consider the writing and reading of a single- 
photon wave packet with a duration T satisfying the con- 
dition Tji <^ T. The writing process will be divided into 
several time intervals, the duration of which is of the or- 
der of Tn and much less than that of the single-photon 
wave packet T. Between the intervals, the atomic system 
is subjected to 27r pulses, which transfer it to different 
subradiant states. The duration of these pulses tp should 
be short in comparison with Tfj. The probability of the 
failure of such a transformation is (see Appendix IB|) 



\2'A-B) = ^i\2A,0B)-\0A,2B)). (21) 



/^failure 



|co| 



5a 



(^^H/2) 



(24) 
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where |coP is the probabihty that the atomic system is 
in the ground state and the field in a single-photon state, 
at the beginning of the transformation, is the Rabi 
frequency of the coherent pulse and ga is the atom-field 
coupling constant on the transition |0j) — Therefore 

if flR ^ (7a, Pfailurc bcCOmCS cloSC tO 0. 

The probability amplitude of photon absorption in 
each time interval is proportional to the wave-function 
amplitude of the photon in accordance with Eq. 
Therefore, at the end of the writing process, the state 
of the atomic system will be determined by a superposi- 
tion of orthogonal subradiant states, each corresponding 
to that the photon was absorbed in the corresponding 
time interval with a weight proportional to the single 
photon wave-function amplitude. If now the atomic sys- 
tem is subjected to a sequence of 2tt pulses, which con- 
vert it back from the subradiant states to the superra- 
diant states, the single-photon wave packet will be re- 
constructed, since the probability amplitude of a sponta- 
neous photon emission is proportional to the amplitude of 
the excited atomic state at the start of the corresponding 
time intervals in accordance with Eq. H12|l . The simple 
example of single-photon writing and reconstruction is 
shown in Fig. |3 In this case the atomic system is di- 



vided into four spatial parts so that the preparation of 
three orthogonal subradiant states is possible, which cor- 
responds to the following Hadamard matrix 



\ 1 -1 -1 1 / 

Let c„ be the probability amplitude that a photon has 
been captured during the n-th time interval. It is conve- 
nient to write the photon state in the following form 

n 

where |1„) is an elementary single-photon wave packet 
propagating through the sample during the n-th time 
interval and /„ is the probability amplitude to find a 
photon in the state |1„). J^n l*-"!^ ~ ^ ^^'^ ideal 
case c„ = /„ (within a phase factor insensitive to n). 
Then the evolution of the atomic-field system during the 
application of 27r-pulses in accordance with the scheme 
(|25|l may be described as follows 



|0)(/l|ll)-^/2|l2) + /3|l3))- 

Cl|lA-B+C-I3)|vac) + |0)(/2|l2) + /3II3)) ^ {ci\1a+B-C-d) + C2|lA-B-C+D))|vac) + |0)/3|l3) ^ 

{ci\1a-B-C+d) + C2\1a+B-C~d) -I- C3|lA-B+C-I5))|vac). 

I 



The first pulse acts on parts BD, the second on BC, and 
the third on BD again [Fig. |2Ia)]. Now the shape of a 
single-photon wave packet is stored as the set of ampli- 
tudes c„. In order to reconstruct the single-photon state 
we apply three 27r-pulses, which act on parts AD, BD, 
and BC, respectively [Fig.Efb)]. It should be noted that 
after the action of the first read-out pulse the phase of 



the emitted wave-packet becomes the same as for the in- 
coming one. Let e„ be the probability amplitude that the 
photon has been emitted between the application of the 
n-th and (n-l-l)-th read-out pulses. Then the evolution of 
the atomic-field system upon the reading of information 
reads 



{ci\1a^B-C+d) + C2\1a+B-C-d) + C3|lA-B+C-D))|vac) 

|0)ei|li) + {c2\1-a+b-c+d) + C3|l_A-B+c+i3))|vac) |0)(ei|li) + e2|l2)) + C3|1_a+s+c-d) |vac) ^ 

|0)(ei|li)+e2|l2)+e3|l3)), 



where e„ — /„ in the ideal case. Moreover, by rearrang- 
ing the coherent 27r pulses at the reading stage, we can 
permute the superradiant states and form a preassigned 
shape of the emitted single-photon wave packet. For ex- 
ample, the application of reading 2tt pulses to the parts 
AC, BC and BD leads to the reversal of a single photon 



wave packet [see Fig.|2Ic)]. Besides, it is possible to mod- 
ulate the phase of the emitted photon at each elementary 
act of read-out by applying 27r-pulses to spatial regions, 
which are complimentary to those used for the writing. 

The number of subradiant states used for the storage 
of a single photon wave packet determines the time reso- 
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FIG. 3: The population of the excited atomic state |c(t)p 
as a function of time t and duration of rectangular incoming 
wave packet Tph- 



Finit) = \/te/tpYi. The term 'quasi-rectangular form' 
means that the leading and trailing fronts of the pulse are 
not shorter than r^. (If this would not be the case the 
one-mode approximation becomes invalid.) Then from 
Eq. ifTT)) we obtain 



FIG. 2: (Color online) Geometry and the temporal order of 
the action of 27r-pulses preparing the subradiant states of an 
extended medium in the case of four spatial parts and three 
possible subradiant states. A filled area indicates the action 
of a 27r-pulse in the corresponding spatial domain: (a) upon 
writing the quantum state of the single-photon wave packet, 
(b) upon read-out without time reversal, and (c) upon read- 
out with time reversal. 



lution of the quantum memory. The reconstructed pulse 
shape in general case is different from the original one, 
but the probability distribution within the time intervals 
between 27r-pulses is the same. Therefore, such a quan- 
tum memory may be useful, for example, for a quantum 
communication using time-bin qubits or, more generally, 
time-bin qudits. The procedure described here for the 
case of a single photon may be extended in a straightfor- 
ward way in order to write and reconstruct a sequence of 
photons. In this case it is necessary to increase the total 
number of orthogonal subradiant states and correspond- 
ing number of spatial parts as was described in Sec. III. 
By rearranging the read-out pulses we can form not only 
a preassigned shape of the emitted single-photon wave 
packets, but also a preassigned sequence of them. 

In order to estimate the efficiency of the writing 
and reading processes consider the solution lfTT|l for the 
case when all atoms are in the ground state at time 
t = and the single-photon wave packet has a quasi- 
rectangular form with a duration Tph and amplitude 








- 1 


exp ^ 











2tr 



, t < Tph, 
t > Tph. 



(26) 



The maximum of \c{t)\ is equal to 0.9 and is attained 
when Tph « 2.5Tfl. The population of the excited atomic 
state |c(i)p as a function of time for the case when the 
incoming wave packet has a rectangular shape is shown 
in Fig. 13 Thus, the efficiency of the writing process for a 
wave packet of duration T ^ tris equal to 81% when the 
27r-pulses separated by 2.5t/j, which is the optimal value. 
Since the duration of the elementary read-out act is also 
equal to 2.5Tfl, the efficiency of read-out of information 
is equal to 1 — exp(— 2.5), i.e., 91%. Therefore the total 
efficiency of the scheme proves to be about 75%. It is also 
possible to find pulse shapes that will enable storage and 
recall with a similar overall efficiency of 75% in a self- 
consistent regime, where the output pulse has the same 
shape as the input pulse. 

However, if each pulse only need to be stored and read 
once the memory can be operated in a nearly 100% effi- 
cient mode. In accordance with Eq. Hll() the probability 
of absorption is maximal for a propagating pulse like 



t tend 

exp — 

TR V 2t7?^ 



t < u 



(27) 
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where iond ^ ^i?- In this case we have 

exp I — , t < tend, 

'WH \ PS) 



exp 



2Tfi 



, t > ^cnd ■ 



Therefore, |c(i)P = 1 at the end of the pulse t = fend- 
Now, if we apply a 27r-pulse to half the atoms at this 
time, a subradiant state is created with unit probabil- 
ity. Upon read-out we obtain (with unit probability) the 
pulse shape 



f ^ ~ ^rcad 



(29) 



corresponding to spontaneous emission after the applica- 
tion of a short read-out pulse at the time t — froad (the 
phase of the emitted photon can be changed by tt at the 
read-out as was considered above). Thus, nearly unit 
efficiency quantum storage with time reversal is possi- 
ble for wave-packets of the form in Eq. 1(23 . Such a 
regime may be useful, for example, for a long-distance 
quantum communication using quantum repeaters |3fl| |. 
when the qubits are only stored and recalled once be- 
fore being measured. Assuming that faint laser pulses 
are used for carrying the information [sj], the prepara- 
tion of a time-bin qubit, which is a superposition of well 
separated exponential wave packets like H27|) . is not a 
challenging experimental task. Each of the wave-packets 
may be recorded (reconstructed) using only one transfor- 
mation to (from) a subradiant state. 

In conclusion of this section, it should be noted that the 
simple model of quantum storage with time reversal may 
be generalized to include propagation effects. In this case 
the shape of pulses to be stored should be determined 
by the time-reversed response function of a sample. As 
a result, both recorded and reconstructed wave-packets 
prove to be modulated in time, but the unit efficiency of 
quantum storage is maintained. 



V. IMPLEMENTING THE SCHEME IN A 
SOLID STATE MEDIUM 

Solids and especially rare-earth-ion-doped crystals are 
attractive materials for optical and quantum optical 
memories. At cryogenic temperatures optical transitions 
of rare-earth-ion-doped crystals have very narrow homo- 
geneous lines, which correspond to long phase relaxation 
of the optical transitions, up to milliseconds (32]. But 
observation of collective spontaneous emission in solids 
as well as related effects such as superradiant forward 
scattering considered above is a nontrivial experimen- 
tal problem. Since inhomogeneous broadening of optical 
transitions is usually several GHz and oscillator strengths 
are small, it is in order to fulfill the condition tjj <C Tj*, 
necessary to use very high densities of impurities and per- 
form observations in the ps regime or to use frequency 



selective excitation of the inhomogeneous profile, which 
significantly limits the superradiant effect. From this 
point the technique of preparing of narrow absorbing 
peaks on a non-absorbing background, i.e. isolated spec- 
tral features corresponding to a group of ions absorbing 
at a specific frequency, in rare-earth-ion-doped crystals 
m HI Hi S |3i il can be very useful. Such specific 
structures can be created as follows. First, spectral pits, 
i.e. wide frequency intervals within the inhomogeneous 
absorption profile that are completely empty of all ab- 
sorption, are created using hole-burning techniques. The 
maximum width of the pits is determined by the hyper- 
fine level separations |37l | . Then narrow peaks of absorp- 
tion can be created by pumping ions absorbing within a 
narrow spectral interval back into the emptied region. It 
is possible to select only a subset of ions absorbing on a 
transition between specific hyperfine levels in the ground 
and excited states '3^. The peaks can have a width of 
the order of the homogeneous linewidth, if a laser with a 
sufhciently narrow linewidth is used for the preparation. 

We now turn to a specific example of how to real- 
ize quantum state storage under the optical subradiance 
regime, using the ^i/4(l) — "'^£'2(1) transition of Pr'^"'' 
ions in Y2Si05. The optical properties of this crystal 
has been studied extensively and the preparation of nar- 
row absorption peaks have already been demonstrated 
on this transition in experiments mentioned above. The 
ground state is split into three doubly degenerate states, 
which are usually labeled [±1/2), |zb3/2), and |±5/2), 
separated by 10.2 and 17.3 MHz |33,|40| for ions in site 
1. A narrow absorption peak with a line- width of less 
than Finh — 100 kHz is usually created in one of these 
states on a 10 — 15 MHz wide pit in the absorption pro- 
file. The life-time of the excited optical state Ti = 164 /is 
and the homogeneous line- width is 2.4 kHz. These val- 
ues correspond to crystals at liquid helium temperature. 
The total inhomogeneous broadening of the optical tran- 
sition is of the order of 4 GHz, so that the Pr ion density 
within a 100 kHz frequency interval is roughly equal to 
2 X 10^" m~^ for crystals with Pr doping concentration 
of 0.05% (9 X 10^** m~'^). The principal point here is the 
possibility of changing this density from that maximum 
value down to zero when preparing the narrow absorp- 
tion peaks. This can be a convenient way to investigate 
a transition from incoherent spontaneous emission to su- 
perradiance in a given crystal by adjusting the super- 
radiant life-time th to the desirable value. Taking the 
wavelength of the optical transition A = 606 nm, the 
length of a sample = 5 mm, the diameter of a focal 
spot 100 /im, we obtain te = 17 ps, fj, = 0.5 x 10~^ 
and Tii = 3.7 ns. The last value is too short because 
the pit width is equal to 10 MHz so that the duration 
of superradiant decay should be longer than 20 ns. If we 
decrease the ion density down to 2 x 10^^ m^'^, we obtain 
obviously = 37 ns, which may be considered as the 
shortest possible value. On the other hand, the upper 
limit for the superradiant decay time is given by the in- 
homogeneous life-time Tj* = (""rinh)""'^, which in this case 
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is equal 3 fis. With the parameters mentioned above we 
obtain Tfi — 1.9 /iS for an ion density of 4 x 10^^ m~'^, 
which may be considered as a lower bound. Thus, we can 
conclude that the observation of optical superradiance in 
the time scale 10^ — 10'^ ns is possible. It should be noted 
that by decreasing the inhomogeneous broadening of the 
absorption peak down to 10 kHz we obtain an inhomo- 
geneous life-time of Tj* — 30 /is. Therefore, in principle, 
it is possible to make tr of the order of several microsec- 
onds, which may be very useful for initial investigations 
of subradiant states in such a system. Since dipole mo- 
ments of optical transitions are rather small (for example, 
0.0078 D for the transition ^Hi{l) - ^02(1) 0), such 
values of allow the use of 27r-pulses with a duration 
up to 100 — 200 ns, thereby reducing the intensity of ex- 
citation pulses. 



A. Active preparation of subradiant states 

The active preparation implies that the phases of 
atomic states are modulated by excitation pulses, such as 
the 27r-pulses considered above. A scheme for preparing 
orthogonal subradiant states that may be simpler from an 
experimental point of view is shown in Fig. 0fa). In this 
case two non-collinear 7r-pulses are applied in sequence 
instead of a single 27r-pulse. Suppose that ki and k2 are 
wave vectors of the first and second pulses, respectively, 
in a pair and the vector ki — k2 is directed along the 
sample z axis. After the excitation the collective atomic 
state |lq) becomes — llq+ki-ka)- It follows that such an 
excitation leads to an added phase shift while propagat- 
ing through the medium. If this phase shift is a multiple 
of 27r, i.e. |ki — k2| — 27rmi^^, where m is an inte- 
ger, a subradiant state is obtained. The rate of collective 
spontaneous emission into the longitudinal mode proves 
to be zero due to destructive interference between vari- 
ous spatial parts. (We propose here that m is small, i.e. 
|ki — k2| ^ loq/c.) Different values of m correspond to 
different orthogonal subradiant states. In order to con- 
vert the system back from the subradiant state into the 
superradiant state the same pulses need to be applied in 
reverse order, i.e. the wave vectors of the first and second 
read-out pulses should be k2 and ki respectively. 

Compared to the excitation by 27r-pulses the advantage 
of such an excitation is that the second 7r-pulse in a pair 
can transfer an atomic state from the excited level to 
another long-lived ground level. In this case the storage 
time may be much longer than the life-time of the excited 
states. 



B. Passive preparation of subradiant states 

On the other hand, it is not necessary to create a phase 
shift using coherent excitation pulses. We can create it 
by controllable phase modulators inserted into the atomic 
ensemble as shown in Fig. ^h). Suppose that each of 



(a) 



Input 



Output 




Writing pair of 7t-pulses 



(b) 



Input 



\b\c\d\ 



Output 



Phase modulators 



FIG. 4: (Color online) Scheme for the preparation of subra- 
diant states in an extended atomic ensemble using (a) pairs 
of non-collinear and non-simultaneous coherent 7r-pulses and 
(b) controllable bi-phase modulators in the case of four spatial 
parts. 



the modulators may be in one of two possible states: an 
off-state, when it does not create a phase shift, and an 
on-state, when a tt phase shift is created. Now, we turn 
on and turn off these bi-phase modulators in a definite 
order instead of application of 27r-pulses. A subradiant 
state is obtained if destructive interference of radiation 
emitted from various spatial parts occurs at the exit of 
the sample and, of course, if the emission process remains 
collective for the whole system. In the case of four spatial 
parts, if we want to obtain the sequence of subradiant 
states which corresponds to the Hadamard matrix H25|l . 
we should switch the phase modulators in the following 
order: ABCD ABCD ABCD, where A {A) means 
that the modulator after the part A is turned on (turned 
off) and so on. The read out without time reversal is 
described by the following scheme: ABCD ABCD 

Abcd. 

Compared to the active preparation scheme such a pas- 
sive preparation scheme may have several advantages. 
Single-photon wave packets propagating in optical fibers 
or wave guides can be stored and reconstructed without 
using transverse coherent excitation. The only condi- 
tion is that photon propagation time through the whole 
system should be shorter than r^j. Further, an atomic 
ensemble (or each optical center) may be placed in a con- 
fined space or artificial structure like a photonic crystal 
such that spontaneous emission into transverse modes is 
inhibited. In this case controlling the emission and cre- 
ating subradiant states for the longitudinal modes may 
enable storage of qubit states in decoherence free sub- 
states in the ensemble during times which may be longer 
than the excited state life-time of a free atom. 
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VI. CONCLUSION 

It is shown that single-photon wave packets with ar- 
bitrary pulse shape resonant with a transition in an en- 
semble of atoms may be recorded and reconstructed in a 
regime of optical subradiance also in the absence of in- 
homogeneous broadening of the resonant transition. Ho- 
mogeneously broadened lines are more prone to super- 
radiance than inhomogeneously broadened lines but the 
technique described here allows the writing and recon- 
struction of time-bin qubits in an optically dense medium 
also when the time-bin qubit duration exceeds the super- 
radiant life-time for the ensemble. This approach could 
therefore open a new window in the search of materials 
for quantum memories. 
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APPENDIX A 

To solve the equations © and © we use the Laplace 
transformation 



y{z)= / e-^'yit)dt, 
Jo 

r-OO 

zy{z)^y{Q)= / e-''y[t)dt 
Jo 



(Al) 
(A2) 



and obtain 

zfk,s{z) - fic,s{0) = -zgk.sV7V(/>(q - k) 

xc(z + «(wo -cj)), (A3) 

zc{z)-c{0) = -i^5k..^^'^*(q-k) 

xfk^s{z-iiuJo~uj)). (A4) 
For the atomic excitation amplitude we find 

5k,.N/^<^*(q-k)/k,.(0) 

- (A5) 



c(0)-*Ek.. 



c(z) 



z — i{u!Q — Uj) 

|gk,.piV|0(q-k)"p 

Z — i{LUo — LU) 



In the Born-Markov approximation we can evaluate the 
sum in the denominator by first converting it to an inte- 
gral 

^ 7^.11-' du dn, (A6) 

k ^ 



and then employing the identity 
1 



lim - 

2-^0 Z — — LO) 



TT (5(a;o — Lij) -\- i- 



P 



with P denoting the principal value. As a result we ob- 
tain 



E 

k.s 



|gk,.piV|0(q-k)|- 

Z — i{LOQ — Oj) 



F + iA, 



(A7) 



where 



is the collective atomic linewidth and A is the collec- 
tive frequency shift, which may be neglected. Here 
fi — 3\^{8t:S)~^ is a geometrical factor |41|, which de- 
scribes the result of the integration 



J df^k ^(d • £k,.)'0'(q - k) = y/id^ 



(A8) 



for a pencil-shaped cylindrical volume with the cross- 
section S, a vector q lying along the axis of the cylinder 
and dipole moments oriented perpendicular to the axis. 
With the help of the inverse Laplace transformation 
yields 



c{t) = c(0) e 



-rt 



5ksVAfr(q-k)/k,.(0) 



F -I- i{ujo — oj) 



(A9) 



It is convenient to define the incoming single-photon wave 
packet shape by the dimensionless photon density in the 
origin of the reference frame 



^ln(t) 



LxLyLz 

~L3 



1/2 



E/k,.(0)< 



(AlO) 



k,s 



Then with the help of approximations (iv-vi) the equa- 
tion (|A9|I may be rewritten as 



c{t) = c(0)e-*/2rfl 
1 



/ drFi„(t-r)e-"/2r^ 
Jo 



Jo 



(All) 



where = (2F)~^ is the superradiant life-time. The 
second term in the brackets obviously equal to zero since 
the interaction between the atomic system and the radi- 
ation field is turned on at the moment t = Q. Thus we 
obtain the solution Eq. (|11|) . 
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In order to obtain the solution for F{t) we integrate 
equation |(SJ), which gives 

/k..(t)-/k,.(O)-»5k,.\/iV0(q-k) / c{t')e-'<^^"-^^''dt'. 

Jo 

(A12) 

Then we substitute this result into the general formula 



LxLyLz 
L3 



1/2 



and take into account that 



dt' e^'-'^o-'^^'' =7r5(wo-w) + i- 



P 



(A13) 



(A14) 



10 loq-uj 
where the imaginary part may be neglected, and 

J dfik ^(d • £k,.)0(q - k) - ^fxd. (A15) 
As a result after some calculations we obtain Eq. p2|l . 

APPENDIX B 

Let us consider a three-level atom interacting with a 
two-mode electromagnetic field. The Hamiltonian of the 
system in the interaction picture is 



H = h{gaRioa + gbR2ib) + H.c. 



(Bl) 



where Rij — \i){j\ is the atomic transition operator (i, 
j — 0,1,2), a, b are the photon annihilation operators for 
modes with frequency wqi and LO12, respectively, ga, gt are 
atomic-field coupling constants. In the photon number 
representation an arbitrary pure state of the system may 
be written as 

2 00 

IV'(0>=X! X! Cj,„„,„Ji)|na,n6), (B2) 

and its evolution in time is given by 

m)) = exp{-zh~^Ht)\^{0)) = C/W|V(0)). (B3) 
The evolution operator U (t) may be written as 

2 



Uit) = J2 4>^J{t)R^. 



(B4) 



where 



ioo = I + gWr-^[cos{Tt) - l]a, 

t>ai = -z.gafl'lT"^ sin(R), 

io2 = 5a5ba^^-2[cos(^^)-l]fo^ 

hi — cos(rt), 

/,i2 = -iV-^ €\ii{Vt)ghh\ 

h2 = I + glbT-^ [cos{rt) - , 



h, = -4 and r = (g^aat +g2^t5)i/2. 
Let the initial state of the system be 

|V(0)) =co|0)|l,a)-^ci|l)|0,a), |coP + |ci|2 = 1, (B5) 



where I a) = exp(-i|ap) X]rb=o ""'("^O is a co- 



,!)-V2| 

herent state with an amplitude a (|ap = (rif,)). This 
means that at the moment t = the atom being in an 
entangled state Co|0)|l, 0) -f- ci|l) |0, 0) with the mode wqi 
starts to interact with a coherent state of the mode W12. 
With the help of equation HB4|I we find 



co{l + gln-^[cos{nt)-l]} 
'icigafl-^ sm{nt) |0)|l,a) 

— icoga^^^ sin(rit) -I- ci cos{flt) |1)|0, a) 
coagagb^~^[cos{D,t) — 1] 
iciagbn~\m{nt) \ \2)\0a). 



where n = [gl + (f]i^/2)2]i/2 and = 2gb\a\ is the 
Rabi frequency of the coherent field. Consequently at 
time tp = Trfl^^ we have 



mp)) = coil-2gln-' 
-ci|l)|0,a) 
~2coagagb^-^\2)\0a 



l,a) 



(B6) 



Thus, the atom proves to be in the upper state |2) with 
the probability 



ga^R 



(B7) 



and returns to the initial state |'(/'(0)) (having a tt phase 
shift for the state |1) amplitude) with the probability 
I— p. For high Rabi frequencies Qn S> ga the probability 
for a successful phase shift operation, l—p, becomes close 
to 1. 
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